Draft version June 9, 2011 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



O 

(N 

oo 
o 

(N 
> 

in 

(N 



o 
o 



X 



LAUNCHING PROTON-DOMINATED JETS FROM ACCRETING KERR BLACK HOLES: 

THE CASE OF M87 

Brezinski^ F. Hujeirat, A.A.^ 

Draft version June 9, 2011 

ABSTRACT 

A general relativistic model for the formation and acceleration of low mass-loaded jets from systems 
containing accreting black holes is presented. The model is based on previous numerical results and 
theoretical studies in the Newtonian regime, but modified to include the effects of space-time curvature 
in the vicinity of the event horizon of a spinning black hole. 

It is argued that the boundary layer between the Keplerian accretion disk and the event horizon is 
best suited for the formation and acceleration of the accretion-powered jets in active galactic nuclei 
and micro-quasars. 

The model presented here is based on matching the solutions of three different regions: i- a weakly 
magnetized Keplerian accretion disk in the outer part, where the transport of angular momentum is 
mediated through the magentorotational instability, ii- a strongly magnetized, advection-dominated 
and turbulent-free boundary layer (BL) between the outer cold accretion disk and the event horizon 
and where the plasma rotates sub-Keplerian and iii- a transition zone (TZ) between the BL and the 
overlying corona, where the electrons and protons are thermally uncoupled, highly dissipative and 
rotate super-Keplerian. 

In the BL, the gravitation-driven dynamical collapse of the plasma increases the strength of the 
poloidal magnetic field (PMF) significantly, subsequently suppressing the generation and dissipation 
of turbulence and turning off the primary source of heating. In this case, the BL appears much fainter 
than standard disk models so as if the disk truncates at a certain radius. The action of the PMF in 
the BL is to initiate torsional Alfwen waves that transport angular momentum from the embedded 
plasma vertically into the TZ, where a significant fraction of the shear-generated toroidal magnetic 
field reconnects, thereby heating the protons up to the virial-temperature. Also, the strong PMF 
forces the electrons to cool rapidly, giving rise therefore to the formation of a gravitationally unbound 
two-temperature proton-dominated outflow. 

Our model predicts the known correlation between the Lorentz-factor and the spin parameter of the 
BH. It also shows that the effective surface of the BL, through which the baryons flow into the TZ, 
shrinks with increasing the spin parameter, implying therefore that low mass-loaded jets most likely 
originate from around Kerr black holes. 

When applying our model to the jet in the elliptical galaxy M87, we find a spin parameter 
a e [0.99,0.998], a transition radius rtr ~ 30 gravitational radii and a fraction of 0.05 — 0.1 of the 
mass accretion rate goes into the TZ, where the plasma speeds up its outward-oriented motion to 
reach a Lorentz factor F G [2.5, 5.0] at rtr- 

Subject headings: Relativity: general. Black hole physics — galaxies: active (M87) — galaxies: indi- 
vidual(M87) — X-rays: binaries (micoqusars) — accretion disks, relativistic jets — 
theory: MHD 



1. INTRODUCTION 

Astrophysical jets have been observed to emanate from 
many astronomical systems such as around young stel- 
lar objects, in binary systems containing compact ob- 
jects as well as in active gala ctic nuclei and quasars (see 
ICamenziiidll2007t lLivioll20Q9l and the references therein). 
Based on astronomical data and theoretical studies, the 
formation of jets in such systems is considered to be in- 
timately connected to the accretion phenomena. These 
data reveal that the very collimated and fast propagat- 
ing jets are found to emanate from systems containing 
black holes with some sort of correlation to their mass 
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(|Fender et al. |[200l . 

In fact a lot of theoretical effort has been made to ex- 
plain the jet-disk or jet-black hole connections, while 
the nature of the interaction between the black hole, 
jet and disk is still not fully understood. A conclusive 
understanding of this interaction would require carry- 
ing out full three-dimensional, general relativistic, time- 
dependent, radiative magnetohydrodynamic calculations 
of dissipative plasmas, using a multi-temperature and 
multi-component plasma description, which is beyond 
the numerical capability of the solvers available to date 
in computational astrophysics. 

Alternatively, based on the combination of numeri- 
cal, observational and theoretical results, we intend to 
construct a theoretical model that would describe the 
BH-jet-disk intera ction properly. A similar app roach has 
been presented by iHujeirat et al.l (p002. 2003.) ; iHujeiratl 




Figure 1. A schematic description of tiie jet-disk-black liole interaction. At large radii the MF in the accretion disk is weak and disordered. 
The MF in the innermost disk region, or equivalently in the boundary layer (BL) between the outer normal accretion disk and the black hole, 
are dragged inwards by the dynamical collapse of the plasma, where they tiecome strong and of large scale topology (blue lines). Angular 
momentum transport in the BL is mediated by magnetic braking, thereby giving rise to the formation of a geometrically thin transition 
zone -TZ between the BL and the overlying tenuous corona. The baryons in the TZ are dissipative and rotate with super-Keplerian velocity. 
The TZ in this model serves also as a runway for accelerating the energetic baryons further by magnetic reconnection of the toroidal flux 
tubes. 



(|2003ll2004l ). though it apphes for the Newtonian regime 
only. The present work is concerned with the extension of 
the above-mentioned Newtonian model into the general 
relativistic regime. Such an extension is necessary, since 
the concerned region of interaction is located in the vicin- 
ity of the event horizon, where the effects of the spin and 
frame-dragging are most prominent. Having performed 
these modifications, we may apply the model to study 
the formation and acceleration of relativistic jets both in 
micro-quasars and in active galaxies. 
This paper is structured as follows: In Sect. [2] we estab- 
lish the underlying assumptions and derive the govern- 
ing equations in Sect. [31 Section [3] is dedicated to the 
derivation of the model and Sect. [5] to a verification of 
its self-consistency. In Sect. [6] we close with a summary 
and conclusions. 

2. BASICS OF JET-DISK-BLACK HOLE INTERACTION 

Consider the innermost part of an accretion disk 
around a central black hole, where the black hole's grav- 
itational force dominates the dynamics of the flow. At 
large distance, the flo w is well described by the stan- 
dard a-disk model (see lShakura fc Sunyaev 1973). Small 
scale magnetic fields (henceforth MF) were found to 
be amplified by the mag neto-rotational instability -MRI 
([Balbus fc Hawlevll 19911 ). which is capable of turning ini- 
tially laminar into turbulent flows. The generated tur- 
bulent eddies have the effect of viscosity, which converts 
a considerable part of the kinetic energy into black body 
radiation that can be observed in the UV and soft X- 
ray bands. Still, turbulent eddies and magnetic recon- 
nection would maintain the magnetic energy to remain 
subequipartitioned with thermal energy, so that the ra- 
tio of the magnetic to gas pressure: j3 = Pmag/Pgas is 
less than unity almost everywhere in the disk. Follow- 



ing iRothstein fc Lovelacd ()1974l ) , weak magnetic fields 
do not hinder accretion, even when they are of large scale 
structure and threading both the corona and the disk. 

However, there is no reason to expect /3 < 1 in the 
innermost part of the disk, where the effect of the spin of 
the central accretin g black hole becomes significant (see 
IPunslv et al. I [20091 and the references therein). MRI, 
Parker instability, reconnection in combination with the 
significant speed up of inflow may lead to the establish- 
ment of a large-scale poloidal MFs, whose corresponding 
energy could be comparable or even l arger than the 
thermal energy of the embedded plasma (jHuieirat et al.l 
I2OOI . In this case, strong MFs in the boundary layer 
(BL) would suppress the generation of turbulence and 
therefore switch off the the source of local heating. Thus 
the BL would appear much fainter than standard disk 
models, or equivalently, the disk appears t o truncate at 
a cer t ain radius close to the central object (jBelloni et al.l 
[2OOOI: [mHeTrat fc Camen"imd1[2000l) . 
Moreover, regions governed by extremely 
strong magnetic fie l ds a re g enerally mat ter-free. 
iMcKinnev fc Gammid (|2004[ ) and[£i:agild (|2008l ) argued 
that such Poynting flux dominated regions may form 
close to the BH and act as a launching mechanism in 
the sense of Blandf o rd fc Znajek extraction process 
(IBlandford fc Znaiekl [19771) . Although ideal MHD 
solvers are incapable of modeling matter-free and mag- 
netic dominated funnels, their simulations show that 
these form at rather high latitudes and therefore most 
likely initially conditioned with no dynamical coupling 
to the accretion flow in the equatorial region. In our case 
however, the collapse-generated strong poloidal magnetic 
flelds in the BL remains confined to the transition zone 
and does not diffuse throughout the corona, so that they 
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Figure 2. The rotational energy of the matter in Keplerian disks Erot is comparable to the gravitational energy, Eg. Part of Erot is 
converted via the MRI into magnetic and turbulent energies i.e. into Emag and Etur, while the rest flows with the accreted matter across 
the transition radius rtr into the BL. rtr is deflned to be the location, where magnetic and thermal energies become comparable i.e., 
Emag ~ Etij. The plasma in the BL undergoes a dynamical collapse, through which the magnetic energy of the poloidal magnetic field 
Ep^'^ increases inwards as r~*. This in turn splits Erot into a major part that goes into the transition zone, E"^^, whereas the rest is 
being accreted by the black hole in the form of kinetic energy Ej^in. Moreover, the strong PMF diminishes Etur and therefore Eth- In the 
transition zone, a part of Ej^ is converted into toroidal magnetic energy EjJ^^^'' which in turn is partially converted into thermal energy 
EJ^'""* through magnetic reconnection to thermal-support the outflowing wind. 



effectively thread a small portion of the corona only. As 
it will be shown in Sec. (4), the time scale for the PMF 
to diffuse throughout the corona is much longer than 
the time scale characterizing the dynamics of the plasma 
in the TZ. Noting that in the the absence of heating 
from below, rotational-u nsupported corona a round BHs 
are thermally unstable (jHuieirat et al.l[2002t ). and that 
thermal conduction across the PMF lines is much weaker 
than along the field lines, we conclude that BH-corona 
most likely are relatively cold voids and therefore are 
thermally and dynamically unimportant for both the jet 
and the disk. 



The question to be addressed here is: what would be 
the most appropriate configuration of plasma flows with 
/3 > 1 in the vicinity of rotating black holes and whether 
such configurations could lead to the formation of the 
highly energetic jets observed to emanate from systems 
containing accreting black holes? 

Let rtr be the radius which separate two regions: ex- 
terior to rtr the flow is said to obey the standard a- 
disk description. Interior to rtr, MFs are predominantly 
poloidal and of large-scale topology, where ideal MHD 
approximation holds. In this case the magnetic field lines 
will be dragged inwards with the plasma particles un- 
der the action of the gravitational force of the central 
BH. As the infiowing plasma rotates, angular momen- 
tum would be extracted and transported vertically on 
the dynamical time scale. This time scale is compara- 
ble to tA = Hd/v^, where Hd is the disk half-thickness 
and v^ is the Alfven velocity along poloidal field lines. 
A larger rtr implies that more rotational energy would 
be extracted from the disk. A too large rtr however, 
would force the inner disk to start collapsing at lager 
radii, hence the inwards-drifted MFs may become suflS- 
ciently strong to terminate accretion completely as in the 



case of " ma gnetically arres ted accretion" , that has been 
reported by llgumenshchevl ([2008). 

Indeed, our model predi ct that rtr < ^9.-'^^' where rs is 
the S chwarzschild radius (jHuieirat et al.l 120031 : iHuieirati 
|2004[) . This is much smaller than the typical dimension 
of the surrounding accretion disc, hence the importance 
of studying the innermost boundary layer (BL) between 
the disk and the central BH. 

Supplied with angular momentum from the BL, the mat- 
ter may start to rotate faster as it moves vertically in 
the manner depicted in Fig. ([3]). On the other hand, 
inspection of the evolution equation of the toroidal mag- 
netic flux implies that the generation of - B^p is propor- 
tional to de^. As this gradient change signs at the in- 
terface betwwen the disk and the corona, toroidal mag- 
netic flux tubes of flux opposite orientation must be gen- 
erated. Part of these tubes will subsequently intersect 
and reconnect, hindering thereby the transport of angu- 
lar momentum to higher latitudes. The trapped rota- 
tional energy in the transition zone (TZ) between the 
disk and the overlying corona forces the plasma to ro- 
tate super-Keplerian and therefore become gravitation- 
ally unbound, hence starts to accelerate outwards. 

The characterizing features of our model compared to 
other theoretical and numerical models are two folds: 
1- Formation of a geometrically thin turbulent-free re- 
gion between the standard disk and the event hori- 
zon, where the plasma components are thermally cou- 
pled but rotate sub-Keplerian, and 2- Formation of two 
geometrically thin runaway zones that sandwich the 
above-mentioned turbulent-free region. The plasmas 
in these zones are proton-dominated, dissipative, two- 
temperature and super-Keplerian rotating. 
We note that simulations that rely on solving the ideal 
MHD equations alone are not not capable of captur- 
ing t hese features properly (e.g. iMcKinnev fc Gamnii^ 
120041) . On the other hand, sophisticated numerical cal- 
culations that relies on using the highly robust 3D axi- 
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symmetic implicit Newtonian rad iative MHD solver - 
IRMHD3 has been carried out by ()Huieirat et al.ll2003D . 
These calculations have explicitly confirmed the forma- 
tion of the above-mentioned features and predicted that 
a fraction of roughly M^/Md ~ 1/20 of the inflowing 
matter goes into gravitationally unbound wind. 

In the present model, the matter in the TZ is pro- 
vided by the inner disk by means of slow vertical drifts, 
v^, mediated by the thermal pressure and MFs. v^ is 
much slower than the other two velocity components. 
Due to the high reconnection rate and inefficient cool- 
ing, the protons in the TZ may easily be heated to reach 
the virial temperature. The super-Keplerian rotating 
plasma in the TZ is centrifugally accelerated outwards. 
A fraction of the toroidal MF in the TZ is advected with 
the wind. At a larger distance the geometrically diverg- 
ing outflow ceases to be diffusive, where Lorentz-forces 
start to redirects the motion of the outflowing proton- 
dominated plasma to form a collimated jet . 
Numerical calculations of ^Huieirat et al.l ()2002| ) have 
shown that a fraction of roughly M„/Md « 1/20 of the 
inflowing matter goes into the form of a wind. 

These properties are the foundation of the present 
model. The subject of the next two Sects, will be to 
deliver an analytic description within the context of gen- 
eral relativity. 




Figure 3. A schematic description of the time-evolution of the an- 
gular velocity Q and its derivative dg^l along the polar coordinate 
8 across the bounday layer for the time sequence ti < t2 < ts < t4. 
The angular momentum of the plasma in the BL is transported pre- 
dominantly in the vertical direction by magnetic braking, thereby 
forming a thin transition zone between the BL and the overly- 
ing tenuous corona, where the angular velocity peaks out. On the 
other hand, since dg^l becomes negative as particles move from the 



disk to higher latitudes, and since B,^ 



gfl, we expect the TZ 



to accommodate a large number of flux tubes of opposite polarity, 
therefore giving rise to an enhanced magnetic reconnection. 



3. THE GOVERNING EQUATIONS 

In this Sect, we quote the equations governing the 
accretion flow. The mass of the central black hole is 
assumed to be much larger than the total mass in the 
surrounding accretion disk and thus self-gravitation of 
the flow is negligible. As a consequence we can describe 
spacetime by means of the Kerr metric. The equations of 
motion of the plasma are derived from the condition of 



conservation of energy, i.e. of the vanishing divergence of 
the stress-energy tensor, V • T = 0, which consists of four 
independent equations. The three spatial components 
yield the momentum equations while the energy equation 
is given by the projection on the timelike four- velocity of 
the plasma, x • (V • T) = 0. No viscous contributions 
to the stress-energy tensor are considered, since, in the 
present model, these are suppressed by strong MFs in 
the inner region of the accretion disk. In the absence of 
viscosity, there are three forms of energy that contribute 
to the stress-energy tensor which can then be written as 

T = Tp + Th + Tem, (1) 

where Tp, Tem and Th represent the stress-energy ten- 
sor of a perfect fluid as well as contributions due to elec- 
tromagnetic fields (henceforth EMF) and thermal pro- 
cesses, respectively. They read as follows: 



(2) 
(3) 



Tem^ - ^ (^F^'F^p ~ \s'^.F''^F^^j . (4) 

c, p and X = dx/dr are the speed of light, the rest- 
mass density and the four-velocity of the plasma with 
proper time t, respectively. T-L = + £ + P/p cor- 
responds to the specific enthalpy, £ to the internal en- 
ergy per mass and P comprises gas, radiation and tur- 
bulent pressure but not magnetic pressure, q is the heat 
flux vector which describes energy fluxes caused by var- 
ious heating and cooling processes. It is purely spatial 
i.e. perpendicular to the fluid worldlines: x • q = 0. 
The most relevant processes in accretion flows are cool- 
ing through bremsstrahlung, comptonization and syn- 
chrotron radiation as well as heating by viscous dissi- 
pation and magnetic diffusivity. Other processes that 
cause heat flux are Coulomb coupling between the elec- 
trons an d protons, adia batic compression and heat con- 
duction ()Huieiratll200l . 

The electromagnetic stress-energy tensor accounts for 
the energy content of the electromagnetic fleld (EMF), 
stresses exerted on the fluid by Lorentz forces and ohmic 
heating caused by electric currents running through a 
resistive plasma. This part is expressed in terms of the 
electromagnetic fleld strength tensor F. In order to take 
EMFs into account the system of equations has to be 
augmented by Maxwell's equations. In terms of F they 
read 



„ 47r 
V • F = j, 

c 

V • *F = 0, 



(5) 
(6) 



where {*F)f^y — ^Sp.upaF'"^ is dual fleld strength tensor, 
j the four-current density and e the totally antisymmet- 
ric Levi-Civita symbol. Before we proceed to the diver- 
gence of the stress-energy tensor, we introduce Ohm's 
law in order to derive the general relativistic induction 
equations. Unfortunately Ohm's law can be very compli- 
cated in_general relativity (for a detailed treatment see 
iMeieii [20041 ). In order to obtain an analytic solution to 
the GRMHD equations, we have to settle for its most 
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simple and straightforward approximation, 



the momentum equations: 



47r , 

— iImJ 
c 



(7) 



where rjM = c^/{4,TTa) is the magnetic difFusivity and a 
the electric conductivity. This is the natural generaliza- 
tion of the Newtonian formula —rjMj = cE + v x B. 
Using Eq. ^ the resistive induction equations may be 
written in the form of an evolution equation for 



dT 



bPF, 



(8) 



where [] indicates antisymmetrisation. Eq. (|8]) contains 
three equations determining the evolution of the electric 
and MF components, respectively. Yet it will be much 
simpler to use Eq. (O to recover the electric fields. 
Introducing the continuity equation 



(px^ = 



(9) 



and making use of eqs. ([S]), © and © one can derive 
the energy equation from i • (V • T) = 0, 

pVif = -pe + —mif + J2 ' (10) 



where = V • a; is the four-dimensional volume expan- 
sion. The first term on the right-hand side of Eq. (jlOp 
accounts for heat generation by compression and the sec- 
ond describes ohmic heating by electric currents running 
through the fluid. The are heating- and cooling func- 
tions that represent all contributions from the heat flux 
vector q, namely 



A-con 

Ape 

^syn 

Ac 



heat conduction 

electron-proton Coulomb coupling 
cooling by synchrotron emission 
Bremsstrahlung-cooling 
cooling by componization 



(11) 



The term ^riMj"^ in Eq. ((TU)) will also be expressed in 
terms of a "heating function, Ao/tm, in regions where we 
the ideal MHD approximation is invalid. 
The three spacelike components of V • T yield the mo- 
mentum equations. In the vicinity of the central object 
gravity is the most dominant force. We will therefore ne- 
glect thermal contributions to the equations of motion, 
i.e. we set H — and g = in the momentum equa- 
ti ons (similar to the condition of negligible speci flc heat 
in lNovikov fc ThorndlTOTllPage fc Thorndll974l l. Using 
Eq. ([9|) we obtain 



c 



(12) 



Now we have derived the complete set of the GRMHD 
equations. In order to perform any speciflc calculations, 
however, we will have to rewrite them into a more 
practical form. Let g be the metric, then the GRMHD 
equations read: 



the continuity equation: 

Videtgl ^ ^ 



(13) 



pX^-df^Xa = ^pxPx''dag,j,u - daP + ^FapjP, 

the energy equation: 

px^d^s = -pe + (±A.) , 

i 

the induction equations: 



(14) 



(15) 



An 



{df, (riMju) - 9^ (?7A/j^)) , (16) 



the electric fleld equations (Ohm's law): 

XtF = -x^F r]Mj ■ 

c 



(17) 



Now the GRMHD equations have been rewritten into a 
form that can be readily applied to any specific model. 
In the following it is assumed that the central object is 
a rotating black hole so spacetime can be described by 
the Kerr-metric. This, however, is not a requirement of 
the model, which is independent of the central object in 
its basic assumptions. Further, we assume that the ac- 
cretion disk lies in the equatorial plane and the system 
is reflection-symmetric. The coordinate system used is 
identical to the pseudo-spherical Boyer-Lyndquist coor- 
dinate system except that the polar coordinate 6 is re- 
deflned by 9 — >■ 7r/2 — 9. In the Newtonian model of 
IHuieiratl ( 20041) the domain of interest was geometrically 
thin, thus we stay close to the equatorial plane: « 0. 
Under these circumstances the Kerr-metric can be ap- 
proximated by 



g = -a^c^dr + {d(p - Ludty + —dr^ + r^dO^, (18) 

where the metric functions in the equatorial plane, cor- 
rect up to order 9^ , read 



VJ 



2r 3a2 



2„2 



2rgr + r,, a 



(19) 



where rg = GM/c^ is the gravitational radius, M is the 
mass of the central object and a G [—1, 1] is the non- 
dimensional Kerr-parameter. 

In the following we use the tetrad system of the zero an- 
gular momentum observer (ZAMO, see e.g. [Camcnzin^ 
I2OO7I ). The components of the EMF as well as the cur- 
rent density, velocity and Lorentz factor are understood 
as measured by ZAMO. For the sake of clarity, however, 
the explicit expressions of these quantities, the sets of 
basis vectors and one-forms are given in the Appendix. 

4. CONSTRUCTING THE COMBINED SOLUTION 

In this Sect, we describe the construction of the com- 
bined solution for both the BL and the TZ in the vicinity 
of the central object. The solution is given based on the 
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time independent and axisymmetric {dt = = 0) equa- 
tions, is taken to be much smaller than v'^,v'^ and 
will be neglected where this is appropriate. 
Let Hd be the disk half thickness in the BL and i/^ the 
thickness of the TZ. Both the BL as well as the TZ are 
assumed to be geometrically thin, i.e. H^/w ^ The 
relative half thickness Hd/r of the BL is further assumed 
to be constant, while the thickness of the TZ is allowed 
to vary with radius, 

Hd/r^smed~9d, (20) 
H^/r = sin 6l^(r) - sin6'd « 9^{r) - dd, (21) 

where 9d is a constant while 6u] is an unspecified function 
of r. Defining the surface densities and E^,, 

Ed = J dOrp^ 2HdP, = j dOrp^ H^p, (22) 

as well as the total mass M inside a given three- 
dimensional, spacelike hypersurface V , averaged over the 
infinitesimal time slice 5t, 

t+st 

M := ^ [ dt [ drdOdipr^ P-, (23) 
St J J a 
t V 

one can derive from the continuity equation. 



^dr (rpVAjv''^ + ^de {piv") - 0, 



(24) 



Md = -27rVAE<j7dwS, 
^ 27r^/AE„7„^;;, 



the well known expressions for the accretion rate, Md, in 
the disk and the outflow rate, in the TZ, 

(25) 
(26) 

where 7 and have been approximated by their vertical 
averages in the integration process. The TZ is fed with 
matter from the disk in the BL by means of a small ver- 
tical drift . Correspondingly Md and M^ are neither 
independent of each other nor constant with radius but 
obey the relation 



drMd{r) = 2drM^{r), 



(27) 



where the factor of , two accounts for the two surfaces of 
the disk. Md and M^ further obey the boundary condi- 
tions 

Md(rt,) = M, Mdire) = Md,B, (28) 

2M„(rtr) = MMtr, M^ire) = 0, (29) 

where M corresponds to the total accretion rate at rtr 
and the relative outflow rate in the TZ is defined by 



M{r) ^2M^{r)/M. 



(30) 



For the moment the _ derivation of Ai and hence the 
profiles of Md and M^ is postponed since it requires 
prior knowledge of several solutions which we have to 
derive first. 



Now we turn to the basic assumptions that iHuieira^ 
made for his Newtonian model. The first assumption is 
that the ideal MHD approximation holds inside the BL 
of the accretion disk. In the non-relativistic limit the 
induction equations read 



^ 47r ^ 
= Vx [v X B TjAij 



(31) 



Correspondingly, the poloidal components are given by 



4:71 



in 



The relativistic generalization of eqs. 
simply given by 



47r , 



(w^v'^B ) do {rjuiv 



An 

-dr {wjv''B ) + —dr {lU-ljip) 



(32) 
(33) 

and ([55)1 is 

(34) 
(35) 



where the covariant component, j^, of the current density 
reads 



47r . 47r , 

—Jv = — K9vvy + 9vt3 ) 



w 



(36) 



In the ideal MHD approximation, the magnetic diffusiv- 
ity vaniushes, ijm = 0, so that 

= de {vo^v^B') , (37) 
= -dr [m^v^'B'^) . (38) 

Eqs. ((34)) . ((35|) allow for a simple solution for S^, 
Bo 



B" 



, where Bq = vj^v^B _ 



(39) 



Thermal equipartition of the poloidal MF at rtr serves as 
boundary condition, thus 



2/zo 



1 ^mp 



Psd(rtr)T,d(rtr), (40) 



where F is the adiabatic index and /i the mean molecular 
weight of the plasma particles. Psd(rtr), Tsdii^tr) corre- 
spond to density and temperature of the standard disk 
at rtr. The radial component B^ can then be derived by 

means of the solenoidal condition ("V • i? = 0") of the 
MF, 

dr {rzuB"-) + de (-B") = 0. (41) 



In the BL, though, the radial component of the MF is 
negligible to the vertical one. Moving to higher latitudes, 
however, the poloidal MF will become increasingly radial 
due to deformation by the wind in the TZ. The profile of 
the toroidal MF can be derived from the toroidal com- 
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ponent of the induction equations, 



^2 



^ ' ^ (42) 



We assume that B'^ is aheady very strong at 9 — 9d so it 
can represent the toroidal field in the TZ while the ideal 
MHD approximation is still valid. Making the approxi- 
mations dr ~ 1/r and de ~ r/Hd, one obtains from Eq. 
dUl) 



LOW 



B-" 



r v:, 



AHdV^d 



-B" 



(43) 

The last term dominates the toroidal field at several 
gravitational radii whereas the first term, containing the 
frame-dragging potential a;, becomes most dominant in 
the immediate vicinity of the event horizon, where the 
poloidal field is predominantly radial. 

In order to construct a reasonable profile for the angu- 
lar velocity, we follow the results of the t h eoreti c al an d 
numerical investigation of IHuieirat et al.l ()2002l I2003D . 
Accordingly, the strong poloidal magnetic fields in the 
boundary layer, i.e. interior to rtr, extract a significant 
fraction of angular momentum from the plasma in the 
equatorial region and forcing it to rotate sub-Keplerian. 
The deposited angular momentum in the TZ enables 
the plasma to rotate super-Keplerian, which susequently 
starts to accelerate its outward-oriented motion on the 
dynamical time scale. In the vicinity of the event hori- 
zon, the dynamical time scale is too short to maintain 
thermal coupling of the electrons with the protons. This 
gives rise to the formation of two-temperature proton- 
dominated gravitationally unbound plasma. The profile 
of the angular velocity in his Newtonian analysis appears 
to fit well to a radial distribution of the form: fl cx r~^/^. 
In the context of general relativity however, the spin of 
the black hole determines the rotational behavior of the 
plasma in the vicinity of the event horizon. Therefore, 
assuming the accretion flow to be in co-rotation with 
the black hole, a reasonable modification of the above- 
mentioned Newtonian profile, would be; 



n 



r5/4rt^.i/4 + r 3/2j 



(44) 



See iBrezins kil (120101) for a detailed discussion of other 
possible profiles. This profile satisfies f2(rtr) = 57k (rtr), 
where riif = V GM ■ (r'^/^ -|-rg^/^a)~^ corresponds to the 
Keplerian angular velocity for co-rotating orbits. The 
modification to counter-rotating orbits is very straight 
forward. 

For fast rotating black holes there will be a radius r* 
where the matter in the BL is non-rotating with respect 
to ZAMO, 



(45) 



where 51 = 51 — a; is the angular velocity as measured by 
ZAMO. Correspondingly, r* is the largest, rational root 
of the equation 



— -2a 



1/4 



0. 



(46) 



Our assumption is that, interior to r.^, the matter in the 
BL is freely falling. 



for r < 



(47) 



which implies that the inflowing matter keeps rotating 
with the frame-dragging frequency w, relative to the co- 
ordinate frame. On the other hand, MFs are still de- 
formed by the frame-dragging effect, thus extracting an- 
gular momentum. The rotational energy is not extracted 
from the matter then, but directly from the central black 
hole (see iPunslv fc Coroniti I fl9m iPunslv et al~l [200l 
for further details). Yet, the total amount extracted in 
this way will be negligible compared to the total rota- 
tional energy of the black hole so no spin-down is taken 
into account. Instead the frame-dragging potential is 
treated as an infinite reservoir of rotational energy. 
In the TZ the matter is rotating with super-Keplerian an- 
gular velocity adopting the same radial profile as in the 
BL. The ri-profiles in the TZ and BL are summarized. 



K 



rtr < r 



(48) 



VGM (r^/Vfii/^ + r3/2a) rs < r (49) 



where the inner boundary tb of the TZ is defined as the 
radius where centrifugal and gravitational acceleration 
are balanced and the effective gravity vanishes, which 
implies that 51 (rg) — 51^^ (r^). 

We will now derive the ^''-profile from the radial 
momentum equation, 



/ r\2 1 2 2/^ (^rP 

-dr 7W ) ^7^1 C C 

Z Z p 



1 

8^ 



dr a^r^B" 



2 Tpr2\ 



dr {AB'^^) dr {r^TU^E 



1 

47rp 



B^deB'^ 



A 

E'-deE'^ 



vjE'^E^'drOJ 



A 



where 



C = -{l-—\dr\na' + dr — 



^V'f^ drfl 



(50) 



(51) 



Regarding that the first term on the right-hand side is of 
order c^/r, one may neglect the pressure term, which is 
of order cs^ /r. Further, we will consider the magnitude 
of the electric field components. From Ohm's law, Eq. 
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(fT7|) . we obtain 

cE'' = v'^B'^ + ^-deB'^, (52) 
7 r 

cE'^ ^v^'Bf ^v'^B'^ ~^^dr(^VAB-^y (53) 

where r/M is the magnetic diffusivity (see Eq. [7]). As 
the plasma motion in the TZ is said to be turbulent and 
dissipative, the finite width of the TZ, Hw, together with 
the Alfwen speed, va, can be used to set an upper limit 
for the turbulent magnetic diffusivity as follows: rju < 
HwVA, where « ^a/i- Inserting rjM in the equations, 
we obtain: 



cE'^ w v^'B'f - yf-B"^ 



(55) 



Since eqs. (f54|) and (|55 
electric fields in Eq, 



show that all contributions from 
(|50l) are smaller than those from 
MPs, it is sufficient to show that all magnetic terms can 
be neglected compared to the net centrifugal and grav- 
itational terms. In this case the centrifugal and gravi- 
tational terms are the only relevant contribution to the 
radial acceleration so they are of order [v^Y /r. Corre- 
spondingly, this assumption is justified if is sufficiently 
larger than typical velocities associated with the electro- 
magnetic contributions. 

It follows from the previous discussion that the MF is 
dominated by B^ in the BL and B'^ in the TZ. Regarding 
that magnetic braking operates on the dynamical time 
scale, i.e. 



(56) 



the magnitude of B^ can be estimated. Eq. ([55]) implies 

that 



= —Vd < ^'d- 



(57) 



Following iGedaliiil ()1993[ ). the relativistic formula for the 
Alfen speed reads 



7A 



(58) 



where 7^ := (1 — w^^/c^) 2. Regarding that % ^ 
and inserting Eq. ([57]) yields 



4:TTPd 



(59) 



Considering that the largest electromagnetic term in Eq 
I 2 



((50|) is of the order Vj^ /r, one can see that all electro- 
magnetic terms can be neglected compared to gravita- 
tional and centrifugal terms within the BL, provided that 

In the TZ the largest electromagnetic term is of the or- 
der v"^ /r. The stationarity condition in the TZ requires 
7^^;^ = 7„v;2 ()Huieiratil200l . Yet, the toroidal field is 
unlikely to be the main driving force since it has turning 
points in its vertical profile. The centrifugal and gravita- 
tional forces on the other hand remain strong throughout 



the TZ. We note that the main effect of EMF both in the 
outer and in the inner part of the accretion disk is to re- 
distribute angular momentum efficiently, so to enable a 
stable accretion of matter. The effect of MFs in the in- 
ner disk is to enable angular momentum exchange of two 
vertically neighboring zones, rather than radially as in 
the case of normal accretion disk. In both cases how- 
ever, MFs have the role of mediator and therefore their 
corresponding total energy must remain well below the 
gravitational and rotational energies. 
Defining the auxiliary function T by 



2v^^ drfl 

c2 -vv^ n ' 



(60) 



(54) one may rewrite Eq. ([50)1 to yield 



drv''=v'"drln a^T 1- 



_(c2_„^2^)^,ln(^a2j-(^l-^jj. (61) 
Eq. ()6ip is a first order ordinary differential equation for 



Its general solution reads 



,,'-2 



7o "o 



J" 



(62) 



where J-' has been normalized to J^j = 1, and where the 
index "0" indicates that the quantity is to be taken at 
radius tq ~ rtr/B in the BL or TZ, respectively. The 
explicit expression of the solution to Eq. (|60|) has been 
moved to appendix |B] It is worth noting that interior to 
r* the matter in the BL is set to be in free fall, implying 
= and J-d = const. The radial velocity then reduces 

to 



r2 
Vd 



a 



J'd.* 1 



""d.tr 



(63) 



We proceed by deriving the profiles of pd from Eq. ([25 



Pd 



Md 



AnHdVA-fdV^^ 



(64) 



where Md — Md{r), and pw from the stationarity condi- 
tion for the toroidal MF in the TZ, namely the equipar- 
tition of kinetic energy with the EMF, ekin,w = cem.w, 
where ekin/EM,w is the corresponding energy density in 
the TZ. Inspection of the 00-component of the stress en- 
ergy tensor ^ yields 



ekin.w = (7^ - 1) PwC^, 



e-EM.i 



Bf"" 
~8^ 



.,r 2 



(65) 



where we have approximated the EMF in the TZ by its 
most dominant components B*^ and E^ and used Eq. 

to set cE'^ « Bfvl,. The density in the TZ then 
reads 

P^ = ^z:2—3, — (66) 
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Ptur/ Pw = Iw^Vtur^, One obtains from Eq. (1571) the rela- 
tion 

(69) 



S. 



For the TZ to be stationary, we require that the amphfi 
cation time scale of B'^, tamp, has to equal the dissipatioi 
time scale tdiss- Inspecting Eq. (|42l) yields 



^amp — 



2 7i« 



ariM 



— tdi. 
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Figure 4. The radial distribution of the velocity of the out- 
flowing plasma in the transition zone at the distance rtr = 40 rg 
from the event horizon is plotted for different values of the spin 
paramter "a". The strong dependence of the wind velocity on the 
spin parameter of the BH in the TZ is most obvious in the vicin- 
ity of the event horizon, where the baryons are strongly accelerated 
outward in order to overcome the deep gravitational well of the BH. 
The dashed line corresponds to the profile of the radial velocity iij 
in the BL, which depends weakly on "a". 

In the TZ the plasma is subject to extensive centrifugal 
forces. We will now investigate what kind of forces are 
opposing this collapse and how the geometrical thickness 
of the TZ can be determined from the vertical momentum 
equation. 



Writing the magnetic difFusivity ryjv/ as 
and using eqs. (|69|) and (|70p yields 



B^ 



rS 



(70) 



(71) 



(72) 



Thus, the thickness of the TZ is directly related to the 
strength of the MF components. Since B^ /B^ <^ 1, we 
conclude that the TZ must be geometrically thin. 



We note that the poloidal magnetic field remains con- 
fined to the TZ and would not diffuse throughout the 
corona. This can be verified by comparing the dynami- 
cal time scale in the TZ (tJ^) to the diffusion time scale 



2 / %'a2v^2 



-B^d^ iraB' 



6*7 
1 

Anp 
where 



1 rga uj 

52 " 



1 



P 



of the PMF into the corona t. 



PMF . 
diff ■ 



r 



J.TZ 

dyn 
-PMF 
^dtff 



{rtr/v. 



Tz) 



mag \ 
'corona ) 



(67) 



(r2-,r,2a2) 



y^2 



2r„aVA 



(68) 

The vertical momentum equation (|67p reveals several 
forces which are due to: 

1. gas pressure Pg 

2. turbulent pressure Ptur 

3. magnetic pressure due to the toroidal MF 

B'^V(87r) 

4. pressure and tension due to the poloidal EMF 

First we neglect v^'^Vg^d? /r^ compared to the cen- 
trifugal term v'^'^ jS"^ in Eq. (1571) . since the former is 
damped by a factor r^^ Super-Keplerian motion ren- 
ders gas pressure as well as all electromagnetic terms ex- 
cept B-^VCStt) negh gible compared to centrifugal forces. 
However, Ptur is most likely responsible for opposing the 
TZ-coUapse. This is due to the fact that the toroidal 
MF has a turning point in the TZ giving rise to re- 
connection of toroidal flux tubes (|Huieira g i200l . Ap- 
proximating de by —r/H^ and 9 « Hu]/r and writing 



where the transition radius rtr bas been taken as a char- 
acteristic length scale both for the dynamical motion 
and for the magnetic diffusion in the TZ, va denotes 
the Alfwen speed. As the centrifugal force acting on the 
super-Keplerian rotating particles in the TZ is the dom- 
inant force and therefore much stronger than the force 
due to magnetic tension that acts to straighten the mag- 
netic force lines, we conclude that tJ^^ <^ '''diYf^ ■ 

Let us now derive the relative outflow rate M from the 
angular momentum equation. 



1 



An pa 



B'-drVA + B'^de] B 



(73) 



In order to do this, we have to integrate over the BL 
in the vertical direction, i.e. from —9^ to 9d- The 
poloidal MF is mainly vertical in the BL and, further, 
at the interface between BL and TZ the vertical rate of 
change of the toroidal MF, dgB"^, exceeds the radial one, 
drB"^ by a factor of order r/Hd so we may savely ne- 
glect B^'dr (yAB'fi^ compared to B'^deB'P. Then, after 
integration, one obtains 



Md drld 



B"Bf 



Using eqs. ([26 



Attw r An 
30)) and combining with eqs. 



(74) 
(I66l) and 
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Figure 5. The two-dimensional distribution of the density, 
poloidal- and angular velocity in the region are displayed, where 
the interaction of the inflow of matter from the accretion disk, the 
accreting black hole and the outflow in the TZ is most effective. The 
poloidal velocity distribution Vp = (v^,v^) in this plot is overplot- 
ted as arrows on the color coded density distribution. The solid 
contours correspond to the angular velocity Q, which attains local 
maxima in the transition layer, TR, where the outward-acceleration 
of the baryons is most significant. 



([72|) for pw and yields 
M = g(r), 



(75) 



g{r) := - aCS 



rc 



where the function C is given in Eq. (|B3p . From this one 
may obtain the expressions 

M.M^M^T^^T' ^W-5M^-4t. (76) 

l-g(r) i-g[r) 

In order to derive the vertical drift through the sur- 
faces of the BL one may act on the continuity Eq. (|24p 

with dO to obtain 



-2rpdidv\ (77) 



where p and 7 have been approximated by their the value 
at the equator, 6* = 0, and — v^{6d) by its value at the 
interface between the disk and the TZ. Approximating 
dr ~ 1/r and using eqs. ([75)1 and (pSj) yields 



(78) 



Following iHuieira^ ()2004[ ) the plasma in the TZ must 
be treated as a two-temperature flow. Due to the low 
density, the time scale for energy exchange is slow com- 
pared to the dynamical time scale. In the disk, on the 
other hand, the density is much higher implying effective 
energy exchange between electrons and protons. There- 
fore we may use a one-temperature description in the 
BL of the disk while we have to consider two seperate 



Figure 6. This plot shows the two-dimensional distribution of the 
MF^ components. The direction of the poloidal field lines is indi- 
cated by normalized arrows along solid lines of constant A^, where 
A = A°'ea is the four-vector potential. Dashed lines mark curves 
of constant |-Bp|, where Bp = (B^,B^). Due to the motion of the 
plasma the magnetic field changes its topology from dipole in the 
outer region into almost monopole in the vicinity of the event hori- 
zon, where the spatial variation of the toroidal field B'^ is most 
significant. Correspondingly, the velocity of the wind in the TZ is 
nearly parallel to the field lines. 

energy equations in the TZ. The energy equation for the 
electrons reads 



A, 



Ar 



and the one for the protons 



Ape - Asyn -Ab-Ac, (79) 



p^v 



^con ~^ -^Ohm ^pe-> 

.(80) 

where the Aj are the heating and cooling functions, intro- 
duced in Sect. [3] The two equations of internal energy, 
eqs. d79l) and (|80l) . are replaced by one single equation, 
where the contribution of ohmic heating, Aohm, may be 
dropped since the plasma is non-resistive in the BL. The 
main heating sources are adiabatic compression and heat 
conduction. The electrons cool mainly by synchrotron 
emission and the protons by Coulomb interaction with 
the electrons. We assume that heat conduction from the 
hot protons in the TZ suffices to compensate for the loss 
of heat. Then the only remaining source of heat is adi- 
abatic compression. Therefore, in the BL, the energy 
equation is given by an ordinary differential equation 
with the simple solution 



Td = Tdivtr) 



A(rtr) rtr 



7dVd 



1-r 



(81) 



where Td and F correspond to disk temperature and adi- 
abatic index, respectively. In the TZ the most domi- 
nant heating process is ohmic heating due to the large 
magnetic diffusivity. Electrons cool effectively by syn- 
chrotron emission. Neglecting other heating and cooling 
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processes yields 

— Aohm ^ ^syn, (82) 

where, based on the previous discussion, Aohm may be 
approximated by 



An 



Fohowing lNovikov fc Thornd (|1973D . A^yn reads 



A...,.=r^('j2;V5:!, (84) 



where T = Saj'^h'^/c, af being the fine-structure 
constant and /e the fraction of unbound electrons to 
baryons. Thus eqs. dM]), ^ and dH]) imply that 



„r2 



r 2 ■ 



(85) 



In order to maintain stationarity, the heating and cooling 
processes in the TZ must operate on the same time scale 
as the advection of B'^, where tadv can be derived from 
Eq. (H^ . The heating time scale is obtained by setting 
theat = p£/ A, where £p = CyTp^^ is the internal energy 
per mass of the protons, Cy being the specific heat per 
mass. This yields 



the 



CvPwTj 



rS _ 

Jp — iadv • 



Aohm OLVi_ 

Correspondingly, the proton temperature reads 



V a 



1) 



+ wj; 



(86) 



(87) 



5. THE FORMATION AND ACCELERATION OF THE JET 

IN M87 

The jet in the elliptical galaxy M 87 is considered to 
emanate from the center, where a 3 • IO^Mq supermas- 
sive black hole is believed to be residing and accreting at 
a rate M = 1.6 • \Q~^MEdd from an optically thin accre- 
tion disk that surroun ds the BH ()Junor fc BirettallTOOl 
iDi Matteo eTaI|[200l . 

The promising and most relevant feature of this jet is 
that the VLBI observations have restricted the size of 
the jet formation region in M87 to r < 70 rg from the 
center of the BH (jJunor fc BirettalllQOSl Tand that a bulk 
Lorentz-factor of the order ^^uik > 3 is found to charac- 
terize the jet-plasma. 

It should be noted, however, that the effect of the stan- 
dard accretion disk in our model is expressed through 
the imposed boundary conditions at the transition ra- 
dius rtr- These conditions apply also to optically thin 
accretions disks, as far as thermal processes are not con- 
sidered. Therefore, without loss of generality, we may 
assume that the disk thickness at rtr does not differ sig- 
nificantly from H « O.lr. 

Inside the BL, i.e. r < rtr, magnetic braking, rather than 
turbulence, is the main mechanism responsible for trans- 
porting angular momentum. This gives rise to the forma- 
tion of a transition zone between the disk and the overly- 
ing corona, where the baryons become gravitationally un- 
bound due to their super-Keplerian rotation, hence start 



to accelerate outwards in the manner shown i n Figs. 
and (|6l ). Additional details can be found in iBrezinskl 



We note that the relative outflow-rate of rest-mass M.q 
may be estimated through the relation: 

M^^M ah^ (88) 

In Fig. ([71 [S]) we show several profile of the Lorentz- 
factor 7n/ and the outflowing rate M.yj oi the plasma 
in the transition zone plotted as functions of the spin 
parameter a and for different transition radii rtr. The 
blue regions contain possible theoretical values that 
surround those revealed by from observations. 
Our model shows that the Lorentz factor of the grav- 
itationally unbound baryons in the transition zone 
correlate with the spin parameter of the giant black 
hole. Similarly, the TZ is best suited for explaining the 
origin and variability of the recently observed TeV and 
7-ray photons from M87, which suggest th at the central 
SMBH must be spi nning at high rate (jWang et al. I 
120081 : iLi et al. 1 12009D . Furthermore, an accretion disk 
that truncate around rtr ~ "^Org and a bulk Lorentz 
factor around 4 appear to be the most probable values 
that fits with observations. Although other truncation 
radii cannot be excluded, we think that rtr ~ 70 rg is 
reasonable in order for the optically thin accretion disk 
to feed the jet with sufficient baryonic matter. 

6. SUMMARY & CONCLUSIONS 

In this paper we have presented a pseudo-analytical 
model for the formation and acceleration of jets in the 
vicinity of rotating black holes under general relativistic 
conditions. The model is based on matching the solutions 
of the time-independent and axi-symmetric GRMHD 
equations in the three regions: the innermost disk, the 
outer standard accretion disk and the transition zone be- 
tween the innermost disk and the overlying corona. 
The main aspects of our model can be summerized as 
follows: 

1. In the outer region, the matter obeys the 
conditions of standard accretion disks. The 
magneto-rotational instability - MRI, is set to be 
the main driver of turbulence, which amplifies the 
magnetic fields on the dynamical time scale, td. td 
itself becomes shorter with decreasing the distance 
from the event horizon. A narrow transition 
region centered around rtr will be established, 
where magnetic and thermal energies become 
quantitatively comparable. 



2. Inside rtr, the collapse of the plasma in the BL 
causes the poloidal magnetic flux to accumulate, 
suppressing therefore the generation of turbulence 
and diminishing heating via turbulent dissipation. 
Such a development has been verified also in the 
3D MHD simulat io ns o f Hgum enshchcv (200j) 
and iPunslv et ahl (|2009f ). Therefore the disk 
ceases to radiate as black body and turns into 
an observationally dark region. In this region, 
angular momentum is transported predominantly 
in the polar direction through magnetic braking 
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Figure 7. The Lorentz-factor — 7uj of the outflowing plamsa in 
the TZ at the distance rtr from the event horizon versus the spin 
of the supermassive BH in M87 is displayed. The profiles show 
the strong correlation of 7u, with the spin parameter of the BH, 
but the moderate correlation with the transition radius rtr- The 
correlation of the Lorentz-factor with the spin parameter manifests 
the argument that jets observed with high Lorentz factor most 
likely emanate from the vicinity of fast rotating BHs. The blue 
region is centered with 7u, and "a" values revealed by observations. 

mediated by torsional Alfven waves - TAWs. 



3. The collapse-induced extraction of rotational 
energy via TAWs in the transition zone,TZ, forces 
the plasma to rotate super-Keplerian, hence 
the baryons start outward-acceleration to reach 
relativistic speeds already at r rtr. The toroidal 
magnetic flux tubes Bt, in the TZ are magnet- 
ically unstable due to their mixed polarity. As 
a consequence, the baryons in the TZ experience 
an enhanced heating and acceleration through the 
magnetic reconnection of these flux tubes. The 
virially hot protons decouple thermally from the 
electrons in the TZ, due to the extensive radio 
emission by the electrons through their gyration 
around the magnetic field. During reconnection, 
both electrons and protons would still emit ex- 
tremely hard photons whilst the bulk of kinetic 
energy will be carried with the protons. 



4. The magnetically-induced collapse of the inner- 
most disk is followed then by an inward drift of 
the transition radius on the local viscous time 
scale. Thus, the plasma in the inner disk starts an 
outside-inside self-illumination through turbulent 
dissipation. 



5. The procedure may repeat itself once the transi- 
tion radius coincides with the event horizon of the 
rotating black hole. 

The parameters of the black holes, namely its mass and 
spin, affect both the location of the transition radius rtr, 
the rate of turbulence-generation in the innermost disk 
and the rate of acceleration of baryons in the transition 
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Figure 8. The integrated flux of the outflowing plamsa in 
the TZ versus the spin parameter of the supermassive BH in M87 
for various transition radii rtr is displayed. The anti-correlation of 
with the spin parameter a indicate that light jets most likely 
originate from around fast rotating BHs, whilst heavy jets from 
slowly rotating or Schwartzschild BHs. 

zone (Fig. 3). The transition radius is found to de- 
crease with increasing the spin "a" of the BH, whilst 
the gamma-factor of the outflowing baryons increases 
with "a". Moreover, the effective surface through which 
baryons enter the TZ shrinks with the spin parameter of 
the BH, implying therefore that light jets must originate 
from around fast spinning black holes, whilst heavy jets 
from slowly rotating or Schwarzschild black holes. 

Furthermore, the plasma in the innermost region is 
expected to enter the phase of turbulent-saturation on 
time scales that anti-correlate with spin of the BH as 
well, thus outlining the role of the spin in enhancing the 
MRI and therefore the generation of turbulence. 

When applying our model to the jet of the elliptical 
galaxy M87, we found that the spin parameter must be 
near unity in order to agree with observations. Conse- 
quently, the central super-massive accreting BH in M87 
must be a maximally rotating Kerr black hole. 
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APPENDIX 

THE GRMHD EQUATIONS 
The set of basis one- forms and vectors Co, used in this work, reads 

e° = acdt, eo = — (9* + udu,) , 

ac 



1 ^ ^ ^« 
e = —=ar, e\ = a, 

VA r 

= rrf^, 62 = -Sg, 

r 

, 1 
e =w {d(fi — Lodt) , es = — 9<y,. 

07 



(Al) 



The components of the field-strength tensor in the coordinate frame may be expressed in terms of the electric and 
magnetic field components E'^, E^, B'^, B^, Bf, measured by the ZAMO, as 

Ftr = -- (E^ + PfdB^) , F*^ = -E\ 

1 



F^e = -ar {E' - PfdB"-) , F*^ = —E", 

^ ' ar 

\/A 

Fre = -y=B^, F^' = ^B^, 

F^r = F^^ = -{B' + PpdE') , 

Fe^ = rc^B^ F'^ = — {B'' - ^fdE^) 

rzu 



(A2) 



where c/Spr) = muj/a is the velocity of a ZAMO with respect to an observer at infinity. Ft^ and hence E'f vanishes 
identically. The components of the electric current density j*, j"^, and j'^, as measured by ZAMO, may be expressed 
as 

Att .( _ adr [rvjE"-) d&E^ 

3 5 ' ' 

c r 

An deBv 



-J 

c r 

47r .g _ drjVABv) 

J ) 
c r 

At, wdrioTB^ dsB^ W^ETdrUJ 

— J — 2 ' ■ K^'^l 

c r rc 

Next we turn to the four-velocity of the plasma, which reads 

X = idt + rdr + Ode + 'fid^ 
= 7 (ceo + v^'ei + w^ej + v'^e^) (A4) 

in the coordinate- and ZAMO frame, respectively. The are the velocities measured by the ZAMO. On the contrary 
to the standard disk solution we have to include a small vertical drift <^v^,v'^. We will neglect wherever possible, 
though. The velocity components and Lorentz-factor are given by 

aVAdt a ' ' dt 

r de dt 1 

V = — — , 7 = a 



adt' ' dr ^1 - v"-'^ /c'^ - v'p'^/c'^' 

(A5) 
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THE PROFILES IN THE BL AND TZ 
We summarize the results obtained in Sect. IH 



^d 

r 



Pd 
Pw 

B'' 

Hw 

Vm 
M 



x/GM (r5/4rt//4 + rg3/2a) ^ ,re[r„rt,.] 
Lu , r < 



VGM (r^/^V V rg3/2a) 



— i^d - ^) 
a 

— (fi^, - oj) 
a 



A 



ip 2 



ip 2 



Td 



9 2 



Id = Jd,tr- 



atr 



2 v 2 

1 c^-t-^^^ 



2 



lw,B- 



2 if 2 

as 1 c - „ 



M, 



B^l + v^lJc' 
87rc2 72-1 

So 



LOW V 
— + - 

avl, V, 



rb 



B- + ^^%B' 



AHdvl 



tt2 

w u. 

—^Iwvl 
rb 

1 - 9tr 



= 9 



1-5 



■Vd9 



Td = Jd(rtr) 



A(rtr) rtr 



i-r 



TUeC^ / 2mpml ^wVw {it - 1) 



T ferS Bv^ 



c2 + vi: 



Cv a > c2 + v--' ■ 



(Bl) 
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The above expressions depend on the six parameters M^, M, a, r, rtr, re- We wiU reformulate them in terms of the 
non-dimensional, scaled variables 

r M . Md 



X = 



and the following constants and auxiliary functions: 

l^ol 



106 G ■ cr. 



g,@ 

Xtr J 



(« 1 — 10 for stellar systems) 



V 



1/2 



id,tr 



= 1? + 10-2 • B'"^2 (i? + 2 a x-^/^ Od Td'^) 



V = l- 



X x^ 



n T -2 -2 / 1 '"K,tr/B 
yd/w = J'd/w Itr/B Oltr/B I J- 2 



(B2) 



■H2 = 10^ • 



To 

Td/w 

Vd/w 



H. 



sd.tr 



rtr 



Od/^ = Xtr/B'^^ + ax-''/'' 



Qd/w = Id/n 



"d/w\ 



{Qd/v.-^'D-'WVd/^-'-l) 



1/2 



l-2aX^rlB^I'x-'l'^°- 

ip 2 



a2 2a2 



W = 1 + 

t9 = e(r-r*) 



1 , r e [r*,rti.] 

, r < r. 



C=-- 



1 



1 _ 



-a^ ^r^ — rr„ 



-+ 



+ ^d 



3H + rg^a^ 



5 = 



2a2 



2a 



pV2(i_vj-V2 + 



+ 1 + -^ 1 + 



(i-v.)-') 



-1/2 



Td/w = exp 



10 



10 10 



i=i " - % 



7 £• t- 



>fc=l 



A; 



+ f - 2a2/o + — ) 



In — ^ In- 
yo - Ui Vi yo 



where y = {r/igY^^ and yi,j/2,2/3 are the three real roots of the equation 



(B3) 
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- yo~'^y^ 



-2y + 2ayo-^ =0, (B4) 



which are located in the interval [—2, 1.5] for all reasonable parameters. Further we have 



2/1 + 2/2 + 2/3 , 

2/4/5 = 7i ± 



2 2/1^ + 2/2^ + 2/3^ (2/1 + 2/2 + 2/3)^ 



± \/ 2/0-2 _ 



2 4 



^"^"fc)'^'"''""^'' ^^^^ 
where n G {6, 7, 8, 9, 10}. In terms of these variables, the solution can be further reformulated as follows 

Oj^^ = 2.0 • lO^s-i • x-^^^Od-^ for r e 
0^2) =4.1. io5.s-i-m-i a a;-3yv-i /or r<n 



fi^ = 2.0- lO^s-i • 0^-1 

=vz/c = x-y^v-y^o^-^T^ 

7^ =25-1/2 g^-l/2-^^-lyi;l/2 

1.2 • 10-3^ • mm-^ x-^ V'^'^ -Ha"^ Qd~^ 
4.4 • 10-^^ • m-2 ^-ig^ ^^-2. 

Qd-^ Td' v„2 w-2 /35-' (1 + • 

•(l-I?g^V„2-^-l) 1 

B"- = 10^ G • TO-i x-i ^6 i?'' T)-^^' Qd^ 

= 106G-m-ia;-iS6Sd-^W-i/2 
B*' = -108G•m-la;-5/4^6^'^2?-^■H2^• 
•Od-^Sd-^rdW-l/2/?5-' 
=10-2rg-x5/4Bv-iD^2O^57S-i^5 
4.4. IQiiHll! .TOa;5/4S¥'-2pg^-i/2^22. 



^1 =3.6- 10-4•m-la;-l/4B62gy-pV2 5^^2-^• 
•Od-' Qd-' Qn, S Td (1 + /3i,2) • 

= = 7 3 . 10-6 . ^-1 a;i/2 ^^2 ^^3 ^^2. 
Te,^ = 7.9 • lO^if • mx5/8 Sg-i S^-i ^^"3/4 7^3- 

■Od 0.-1/2 5-1/2 7-^-1 ^^1/2 -^^-3/2 y^;5/4. 

(1 - vg^v^'w-'Y^' (1 - /j:;,^)'/' (i + pi^'y 

Tp^y, =7.3 -1012/4:. D-2G^-3/2v^-3>v2. 

• (1 - -Dg^vjw-^) (1 - /3;2) (1 + 2)-i 



